PROCEEDINGS 

0F SCIENCE 



Non-perturbative quark mass dependence in the 
heavy-light sector of two-flavour QCD 



~%LPHA 



CERN-PH-TH/2008-207 
MS-TP-08-25 

Collaboration DESY 08-148 

00 SFB/CPP-08-85 

o 
o 

<N Michele Delia Morte 

CERN, Physics Department, TH Unit, CH-1211 Geneva 23, Switzerland 
E-mail: dellamor@mail . cern . ch 

^ Patrick Fritzsch* \ Jochen Heitger 

Westfalische Wilhelms-Universitat Miinster, Institutfiir Theoretische Physik, 
Wilhelm-Klemm-Str. 9, D-48149 Miinster, Germany 

E-mail: f ritzsch@uni-muenster . de, heitgerguni-muenster . de 

o | Rainer Sommer 

Deutsches Elektronen- Synchrotron DESY, Zeuthen, 
t Platanenallee 6, D-15738 Zeuthen, Germany 

E-mail: rainer.sommer@desy.de 

> 

We present preliminary results of the non-perturbative heavy quark mass dependence of heavy- 
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1. Introduction 



The Heavy Quark Effective Theory (HQET) Lagrangian, 



^hqetO) = VhO) Aj + Wl-fflkinD^OfcpinffB i(/ h (i) + O (l/m 2 ) 



(1.1) 



provides an expansion of QCD amplitudes in the inverse heavy quark mass, 1 / m, and is renormal- 
izable at any finite order in l/m by means of power counting. It is a standard phenomenological 
tool which simplifies the QCD dynamics in the limit of large masses like that of the c- or b-quark. 
The HQET Lagrangian as written in (1.1) consists of a leading, static term, Do, which describes 
the dynamics in the limit of infinitely heavy quark mass, and of the two subleading kinetic and spin 
terms whose coefficients are o\i n = ft) sp j n = 1/ (2m) in the classical theory. 

To explicitly make HQET an effective theory of QCD requires matching calculations to ex- 
press the parameters in the Lagrangian (m, (Qkm, ■ • • ) by those of QCD. In principle one could 
employ perturbation theory in the matching step, but owing to the difficulty of a reliable error 
estimation it may be hard to disentangle and quantify deviations coming from higher orders or 
non-perturbative effects in particular, since there are power-divergent mixings [1], Despite its phe- 
nomenological success — reflecting, e.g., in determinations of V c b [2, 3] or HQET hadronic matrix 
elements where perturbative HQET enters — , independent and non-perturbative tests of HQET 
may provide a deeper insight into the feasibility of the effective theory approach. In the quenched 
approximation, such a test was performed and discussed in [4]. 

In the course of the non-perturbative matching between QCD and HQET, originally proposed 
in [5] and currently being applied to QCD with two (massless) dynamical sea quarks [6] (see [7] 
for the variant of the HMC algorithm we use), we here give a preliminary report on an extension 
of the quenched study [4] to the two-flavour case. The non-perturbative large-mass behaviour of 
some meson observables, computed in the continuum limit of finite-volume QCD, is confronted 
with the static theory in order to investigate the range of validity of the 1 / m-expansion and to 
estimate the size of the 1 /m-corrections. The use of & finite volume of about (0.5 fm) 4 is crucial in 
this context, since clean non-perturbative comparisons of QCD and HQET in the continuum limit 
require m <C 1 /a prior to a — ► 0. 

2. Observables and their large-mass asymptotics 

Our observables are built from relativistic heavy-light Schrodinger functional (SF) correlation 
functions, /a is a correlator between a heavy-light pseudoscalar boundary source and an axial 
current operator insertion in the bulk, f\ a boundary-to-boundary correlation function and k\,k\ are 
their vector channel analogues. We further use the pseudoscalar, /p, and tensor, kj, correlators to 
improve the bare currents in the respective channels. More details and any unexplained notation are 
found in [4]. From now on, these correlators are referred to as their 0(a) improved lattice versions, 
which amounts to replace the quark bilinears such as the heavy-light axial current according to 



->Z A [1 + 2^A(«w qi i + am q) h)] xA M . 



(2.1) 



2 



NP quark mass dependence in the heavy-light sector of two-flavour QCD 



Patrick Fritzsch 



We take the non-perturbative values of Za,Zv from [8, 9] and the 1-loop perturbative estimates 
of [10] for bx, X G {A, V}. Our test observables built from these renormalized correlators are 

T av (L,M) = I [r PS (L,M) + 3r v (L,M)] , R, P UL,M) = \Hf,/k Y ) , (2.2) 

F PS (L,M). + ^, F V (L,M) ee -M^), (2.3) 



where the following definitions of the pseudoscalar and vector effective energies apply: 

fk(T/2) 



r PS (L,M)== -/-ln[/ A (jc )] 
dxo 

r v (L,M) = -— ln[£ v (*o)] 
dx 



. X0= r/2 /a{T/2) 
K(T/2) 

x =T/2 ky(T/2) ' 



(2.5) 
(2.6) 



In physically large volume, as L — > oo, r av becomes proportional to the spin-averaged mass of 
the heavy-light meson, /? S pin to the spin-splitting term and }ps (Yy) to the pseudoscalar (vector) 
heavy-light meson decay constant; in this sense we will also use the shorthand decay constant. 

These quantities are expected to be described by a power series in l/z with logarithmic modi- 
fications, where 

z = L l M, M= lim {[2b g 2 (pL)]- d °l^m(pL)\ , (2.7) 

N f = 2, b = {n-lN i )/(4n) 2 , d = S/(4n) 2 , 

and M denotes the renormalization group invariant (RGI) mass of the heavy quark flavour. Our 
choice of lattices with spatial extent L = L\,T = L, and the corresponding simulation parameters 
are summarized in table 1. In physical units, L\ is about 0.5 fm. For a more detailed account on 
how the z-values under investigation, z £ {4,6,7,9, 11, 13, 15, 18,21}, can be kept fixed for given 
L\/a in dependence of the bare parameters, the reader is referred to [6]. 

As we are interested in the large z-asymptotics of our observables, one also needs their counter- 
parts computed in the effective theory, as long as the latter are non-trivial. We denote the associated 
quantities by a superscript 'stat', e.g. /a —* f^ 1 - I n case °f tne decay constant we then construct 

so that 

limFp S (L,M) =X(L) = lim Yy(L,M) (2.9) 

7 ^OO 7 }CO 

holds at the classical level. Due to the heavy-quark spin symmetry in the static limit, Yp§ and Yy 
converge to the same limit. 

In quantum theory, the scale dependent renormalization of the effective theory introduces the 
mentioned logarithmic modifications. An example is the axial current in the effective theory, where 
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Li /a 


J8 


Kj W Kc 


L\m\ 


20 


6.1906 


0.135997290 


+0.00055(13) 


24 


6.3158 


0.135772110 


-0.000145(66) 


32 


6.5113 


0.135421494 


+0.000143(36) 


40 


6.6380 


0.135192285 


+0.000024(24) 



Table 1: Simulation parameters that correspond to a line of constant physics characterized by g 2 (L\ ) = 4.484 
and L\m\ ss in the light (i.e. sea) quark sector. Suitable choices for the hopping parameter of the heavy 
quark, K"h, allow to fix its dimensionless RGI mass, z = L\M, to a set of desired values extending from the 
charm to the bottom quark region. See [6] for details. 



Z^ at (/x)Xb are (Z,) depends logarithmically on the chosen renormaliza- 



the renormalized Xr(L,jU) 
tion scale jj. as well as on the renormalization scheme. As for the mass, eq. (2.7), this dependence 
is removed explicitly by passing to the RGI matrix element 



Xrgi(L)= Urn {[26og 2 (M)r W(2feo; 



ZRGI^bare(^) . 



7o 



l/(47i 2 ). (2.10) 



For TYf = 2, the renormalization factor is known non-perturbatively from [11]. After express- 
ing our QCD test observables through the corresponding RGIs, their large-z behaviour is driven by 
the RGIs of the effective theory together with so-called conversion functions, C, which contain the 
full logarithmic mass dependence of the associated operators. As arguments of the latter we choose 
the ratio of RGIs M/A, since it can be fixed on the lattice without perturbative uncertainties [12]. 
The 1 /z-expansions of our test observables now read as follows: 

u- 



Y X (L,M) 
Rx(L,M) 

flspm(i,M) 

LT av {L,M) 



M- 



M- 



M- 



Cx{M/A m 
Cx(M/A m 



1 



i+o(i/ z ; 



spin 



yspin 
RGI 



l+0(l/z) 
L)/z 



X = PS,V, 
X = PS/V,PS/P, 



l + 0(l/z) 



M/A 



MS 



l + 0(l/z) 



(2.11) 
(2.12) 
(2.13) 
(2.14) 



Beside the overall logarithmic mass dependence governed by the conversion functions we enclose 
the leading matrix elements in the effective theory in square brackets. The numerical evaluation of 
the Cx proceeds as explained in [4, 11], where always the highest available perturbative approxima- 
tion of the involved RG functions is employed, in particular the 3-loop anomalous dimensions y of 
the axial current and the chromomagnetic operator in HQET, respectively [13, 14]. By comparing 
different loop orders for their evaluation, we conclude that the remaining perturbative uncertainty 
of the conversion functions originating from unknown higher orders is much smaller than the pre- 
cision of our lattice data so that a study of the 1 jz n -corrections becomes feasible. 

3. Results 



As examples for the observables under study, continuum extrapolations linear in (a/L) 2 of 
Lr av and /? sp in for the phase parameter 1 6 = 0.5 are presented in fig. 1. In our volume of extent 

'The SF boundary conditions in space on the fermion fields are taken to be periodic up to a global phase 6. 
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Figure 1: Left: Continuum extrapolations linear in (a/L) 2 of the spin-averaged mass Lr av (z, — 0.5) from 
z = 4 (bottom) to z = 21 (top). The error of the continuum limits is smaller than the symbol size. Right: 
Continuum extrapolations linear in (a/L) 2 of the spin-splitting R sp m(z, 9 = 0.5) from z = 4 (top) to z = 21 
(bottom). — Depending on the value of z, the coarsest or the coarsest two lattices are omitted from the fits. 

L = L\ = T m 0.5 fm, which admits to reach heavy quark masses up to M 1.3Mb, the a — > 
extrapolations appear to be well controllable provided that one accounts for the growing (heavy) 
quark mass in lattice units at given a/L as z is increased. Similar to the quenched work [4], this is 
realized by imposing a cut on aM {aM < 0.7) that translates, for given z, into the coarsest resolu- 
tions which may still be included in the continuum extrapolations. For all observables considered, 
the growing deviation, as z is increased, of the extrapolating fit function from the result at the 
respective coarsest lattice resolution available suggests that between z = 18 and z = 21 the 0(a) 
improvement and thus the ^-expansion has broken down for our lattices. 

Polynomial fits in 1 fz of the continuum Lr av / (zC mass ), (zR sp m) /C S pin as well as of the finite- 
volume pseudoscalar and vector decay constants and their ratio R^s/w (divided by the correspond- 
ing conversion functions) are shown in figs. 2 and 3. If the conversion functions Cx are evaluated 
including the highest available perturbative accuracy, the z-dependence of r av and R PS / v is consis- 
tent with the leading term in the 1 /z-expansion, which in case of /?ps/v i s fixed by the heavy-quark 
spin symmetry of the static theory. Note in addition that the l/z—>0 limit of these observables is 
independent of the choice for the periodicity angle 8 (= 0,0.5, 1 here), as it should be. 

In general, the 1/z-corrections of our observables are reasonably small. As can be inferred 
from the l.h.s. of fig. 3, Fps/Cps and Fy/Cy converge to the same limit (see (2.9)) as expected. In 
principle, this common limit of their z-dependence can be constrained by the static theory as well, 
because the renormalization factor of the associated static axial current matrix element, (2.10), is 
already known non-perturbatively [11]. The computation of the bare matrix element to obtain the 
static result for Xrgi is in progress. 

4. Conclusions 

We have presented a status report of our ongoing study of the large-mass asymptotics of non- 
perturbatively renormalized heavy-light meson observables in finite-volume two-flavour QCD. 
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Figure 2: Left: Asymptotics of ir av /(zC ma ss) versus l/z. All data sets are extrapolated to the (known) 
static limit (= 1, cf. (2.14)) with a constrained quadratic fit. Right: 1/z-dependence of (z^ S pin)/C S pin with 
an unconstrained quadratic fit to all data points. — Blue, red and green symbols refer to 9 € {0, 0.5, 1}. 
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Figure 3: Left: Asymptotics of Yx/Cx versus l/z and its unconstrained quadratic fits to the static limit. Data 
sets approaching l/z = from above (below) corresponds to X = PS (X = V). Right: ^ps/v/Cps/v versus 
1 jz and constrained quadratic fits to the static limit (= 1, cf. (2.12)). 



Their behaviour in dependence of the inverse RGI heavy quark mass is well compatible with the 
predictions of HQET confirming the feasibility of a precise non-perturbative matching of QCD 
and HQET. The only perturbative uncertainties owing to the leading logarithmic mass dependence, 
which is induced by the conversion functions C relating our observables to the RGIs of the ef- 
fective theory, are under reasonable control. It appears that for the studied observables the power 
corrections dominate over the perturbative ones in the considered range of z. 

Our continuum extrapolations may be further improved by a removal of perturbative cutoff 
effects prior to the extrapolations. 
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